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Abstract. A nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used 
to describe the properties of Tycho's SNR. Observations of the expansion characteristics and of the nonthermal 
radio and X-ray emission spectra, assumed to be of synchrotron origin, are used to constrain the overall dynamical 
evolution and the particle acceleration parameters of the system, in addition to what is known from independent 
estimates of the distance and thermal X-ray observations. It is shown that a very efficient production of nuclear 
cosmic rays, leading to strong shock modification, and a large downstream magnetic field strength Bd ~ 240 /iG 
are required to reproduce the observed synchrotron emission from radio to X-ray frequencies. This field strength 
is still well within the upper bound for the effective magnetic field, consistent with the acceleration process. 
The 7r°-decay 7-ray flux turns out to be somewhat greater than the inverse Compton (IC) flux off the Cosmic 
Microwave Background at energies below 1 TeV, dominating it strongly at 10 TeV. The predicted TeV 7-ray 
flux is consistent with but close to the very low upper limit recently obtained by HEGRA. A future detection at 
e 7 ~ 10 TeV would clearly indicate hadronic emission. 

Key words. Acceleration of particles - Radiation mechanisms: non-thermal - Stars: supernovae: individual: Tycho's 
SNR - Radio continuum: ISM - X-rays: ISM - Gamma rays: theory 



1. Introduction 

The observations of Tycho's supernova remnant (SNR) 
with the HEGRA stereoscopic system of imaging atmo- 
spheric Cherenkov telescopes (IACTs) on La Palma have 
recently been analyzed. The object has long been consid- 
ered as a prototype candidate hadronic CR source in the 
Northern Hemisphere (e.g. Drury et al. 1994 ), although it 
was always clear that the sensitivity of the present gener- 
ation of IACTs is marginal for a detection. In fact, after 
~ 65 hours of observation time, HEGRA did not detect 
Tycho's SNR, but it could establish a very low 3cr upper 
flux limit of 5.78 x 10 -13 photons cm~ 2 s _1 , or 33 milli- 



Crab, at energies > 1 TeV (Aharonian et al. 2001). This 
value is about a factor 4 lower than the one previously 
published by the Whipple collaboration (Buckley et al. 
1998] ), assuming a spectral index of —1.1 for the compar- 
ison. In the above HEGRA paper on Tycho the existing 
radio and X-ray synchrotron observations were used to in- 
fer a lower limit to the mean magnetic field strength in the 
remnant due to the nondetection of Inverse Compton (IC) 



emission off the Cosmic Microwave Background (CMB). 
At the same time, published estimates of the hadronic ir°- 



decay gamma-ray emission (Berezhko & Volk 1997) were 
employed and scaled to the parameters of Tycho's SNR, to 
compare with the new upper flux limit. The predictions 
from the time-dependent kinetic model (Berezhko et al. 



1996) were also renormalized to take the expected devi- 
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ations from spherical symmetry of the nucleon injection 
rate for a Type la SNR into account, and were found to 
be consistent with the present HEGRA nondetection in 7- 
rays, although the predicted flux values for the hadronic 
emission were only slightly smaller than the observational 
upper limit. 

This tantalizing situation has prompted us to model 
the acceleration of both electrons and protons in detail 
with the nonlinear kinetic theory, using the observed syn- 
chrotron emission as a constraint on the electron acceler- 
ation characteristics, and thereby to model the hadronic 
and IC 7-ray emission simultaneously. 

We demonstrate that, together with the renormaliza- 
tion, the existing data are consistent with very efficient 
acceleration of CR nuclei at the SN shock wave which con- 
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verts a significant fraction of the initial SNR energy con- 
tent into CR energy. This energy is distributed between 
energetic protons and electrons in a proportion similar to 
that of the Galactic CRs. Therefore Tycho's SNR might 
indeed be a typical Galactic CR source. Direct evidence for 
the inferred strong production of nucleonic CRs in Tycho 
would not only require a 7-ray detection as such, but in 
particular a spectrum measurement, ideally beginning at 
photon energies of about 100 MeV up to the highest en- 
ergies. Of particular importance would be measurements 
in excess of 10 TeV particle energy which are beyond the 
reach of the accelerated electron component that suffers a 
cutoff at energy of about 50 TeV due to synchrotron losses 
in our model. This specific spectral behavior implies a high 
magnetic field strength of about 40 /iG upstream of the 
outer SNR shock. Such a field strength is indeed required 
to fit the nonthermal radio and X-ray synchrotron emis- 
sion. We shall give an upper bound on the effective mag- 
netic field in the discussion section. It well encompasses 
the deduced field strength. 

2. Model 



A supernova (SN) explosion ejects a shell of matter with 
total energy E sn and mass M e - y During an initial period 
the shell material has a broad distribution in velocity v. 
The fastest part of these ejecta is described by a power law 
dM e j/dv oc v 2 ~ k (e.g. Jones et al. |198l| , Chevalier |1982Q . 
The interaction of the ejecta with the interstellar medium 
(ISM) there creates a strong shock which accelerates par- 
ticles. 

Our nonlinear model (Berezhko et al. 



1996; Berezhko 



& Volk 1997) is based on a fully time-dependent solu- 
tion of the CR transport equations together with the gas 
dynamic equations in spherical symmetry. Regarding the 
coefficients in the transport equations and the boundary 
conditions as well as the approximation of spherical sym- 
metry, we shall introduce a number of approxima tions as 
follows (see also the discussion in Berezhko et al. 2002 ): 
The CR diffusion coefficient is taken as the Bohm limit 



k(p) = K(mc)(p/mc) J 



(1) 



where k(toc) = mc 2 / (3eB), e and m are the particle 
charge and mass, p denotes the particle momentum, B 
is the magnetic field strength, and c is the speed of light. 
This limiting value is consistent with the high, probably 
turbulently amplified magnetic field strength (Lucek & 
Bell 2000) due to the strong pressure gradient of the ener- 
getic particles in the very strong shock of the young rem- 
nant. In fact, Eq. (1) uses the relativistic form of Bohm 
diffusion, equal for both electrons and nuclei. In fact, the 
actual form of n(jp) at low momenta p <C p m ax does not 
influence the particle spectrum if n(p) is a strongly in- 
creasing function of p. In addition at nonrelativistic proton 
energies p < m p c, even if the electron and proton diffusion 
coefficients are different at the same p, their spectra f(p) 
are essentially the same because these low energy parti- 
cles do not produce any shock modification and therefore 



their spectrum is f(p) oc p~ q % where q s is the spectral 
index appropriate for the subshock compression ratio a s , 
calculated below. 

The number of suprathermal protons injected into the 
acceleration process is described by a dimensionless in- 
jection parameter r\ which is a fixed fraction of the ISM 
particles entering the shock front. For simplicity it is as- 
sumed that the injected particles have a velocity four times 
higher than the postshock sound speed. It is expected that 
ion injection is quite efficient at the quasi-parallel por- 
tion of the shock surface, where it is characterized by the 
value rj ~ 10~ 4 (see Volk et al. |2002| for details). In mo- 
mentum space the accelerated (nonthermal) proton dis- 
tribution function grows smoothly out of the downstream 
thermal distribution (e.g. Malkov & Volk |l995| ). 

Assuming nucleon injection to occur through the leak- 
age of suprathermal particles into the upstream region 
from behind the shock, the relevant injection velocity par- 
allel to the magnetic field is 14cos -1 </>, where V s denotes 
the (radial) shock velocity and 4> is the angle between the 
shock normal and the downstream magnetic field direc- 
tion; it is refracted away from the shock normal. Therefore 
the injection is expected to be strongly suppressed at the 
quasi-perpendicular surface fraction of the shock and this 
surface fraction is larger than that of the quasi-parallel 
fraction due to the field refraction. This lack of symmetry 
in the actual SNR can be approximately taken into ac- 
count by a renormalization factor / re < 1/2 which dimin- 
ishes the nucleonic CR production efficiency as calculated 
in the spherical model, and all effects associated with it. A 
more detailed estimate, taking also the diffusive broaden- 
ing of the quasi-parallel area into account, yields roughly 
/ re = 0.15 to 0.25 (Volk et al. p002| ). 

This is not necessarily in conflic t with the radio po- 
larization result of Reynoso et al. ( 1997 ) and Dickel et 
al. ( 1991 ) who obtain from their VLA measurements a 
roughly radial ordered magnetic field component over 
most of the shock surface. The fractional polar ization of 
10% to 15% is rather small (Dickel et al. |l99l| ), indicat- 
ing a rather small excess of ordered radial magnet ic fiel d 
(see also the discussion in Reynolds & Gilmore 1993 ). 
Obviously Tycho's SNR is not a magnetic monopole. To 
resolve this puzzle we argue that there exist at least two 
physical processes which tend to produce radial field di- 
rections in a remnant that is still in the sweep-up phase 
and as a whole embeded in a uniform ISM and mag- 
netic field. First of all, such young SNRs as Tycho's are 
Rayleigh- Taylor unstable during the sweep-up phase (Gull 



1973; Velazquez et al. 1998). Secondly, the ejecta them- 



selves may be so nonuniform that they can wrinkle the 



blast wave surface (Aschenbach et al. 1995). In addition, 
the ambient ISM around Tycho appears to be unexpect- 
edly nonuniform as well, cf. Reynoso et al. (1997). Even 



though these irregularities are of a macroscopic, hydro- 
magnetic nature, their size spectrum probably extends 
to quite small scales. The acceleration scales can even 
be much smaller. In fact, the acceleration scales of pro- 
tons with multi-GeV energies, on whose strong self-excited 
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magnetic irregularities (Bell ^978) the electrons may para- 
sitically accelerate as well, is of order 10 -4 times the shock 
radius or even smaller, and the proton scales are connected 
with quasi-parallel, i.e. quasi-radial field line portions. On 
other portions no acceleration occurs, and they will not 
be illuminated by synchrotron emitting electrons. Thus it 
is possible that the magnetic field appears essentially ra- 
dial over the remnant surface even to the VLA. We leave 
it open here to which extent small-scale magnetic irreg- 
ularities that lead to measurable polarization effects can 
also be nonlinearly produced by the accelerated particle 
population itself. 

Due to their small mass m e , suprathermal electrons 
cannot resonantly interact with the hydromagnetic wave 
turbulence created by the energetically dominant accel- 
erating ions. Therefore electrons can only participate in 
the diffusive shock acceleration process at high energies 
< m c c 2 , albeit still nonrelativistic as far as nuclei are 
concerned. Presumably they are inj ected by electrostatic 
fluctuations (e.g. Malkov & Drury 2001 ). For simplicity 
we can consider their acceleration to start from the same 
momentum as protons start, if we restrict our further con- 
sideration to relativistic electrons, at energies possibly far 
above their actual injection energy. What counts alone for 
the considerations in this paper is that at energies corre- 
sponding to those of relativistic protons - only those en- 
ergies play a role for the observed synchrotron emission - 
electrons have exactly the same dynamics as the protons. 
Therefore, below the momenta where synchrotron losses 
become important, the electron distribution function has 
at any given time the form 



f e (p) = K ep f(p) 



(2) 



with a factor K ep that is of the order of 10 



2 for the 
average CRs in the Galaxy. 

Clearly, from the point of view of injec- 
tion/acceleration theory, we must treat K cp together 
with B as free parameters, and r\ as well as / rc as 
theoretically not yet accurately calculable parameters, 
to be quantitatively determined by comparison with 
observations. 

The electron distribution function f c (p) deviates only 
at sufficiently large momenta from Eq. (2) due to syn- 
chrotron losses, which are taken into account by supple- 
menting the ordinary diffusive transport equation by a ra- 
diative loss term. The solution of the dynamic equations 
at each instant of time yields the CR spectrum and the 
spatial distributions of CRs and thermal gas. This allows 
the calculation of the expected flux F^(e 7 ) of 7-rays from 
7r°-decay due to p — p collisions of CRs with the gas nuclei. 
The choice of K ep allows one then to determine the elec- 
tron distribution function in the energy region where losses 
can be neglected, and to calculate the associated emission. 
Details about this calculation are given in Berezhko et al. 



3. Results and Discussion 

Tycho was a type la SN. Therefore we use typical SN 
la parameters in our calculations: ejected mass M e j = 
IAMq, k = 7, and a uniform ambient ISM with hydro- 
gen number density N-r = 0.5 cm~ 3 and temperature 
T = 10 4 K. Following Chevalier et al ( |l980| ) (see also 
Heav ens 1984 ; Reynolds & Ellison 1992 ; Aharonian et al. 
|200lD we in addition adopt a distance d — 2.3 kpc, and 
a present radius of ~ 4'. These parameters are similar to 



those inferred by Dickel & Jones (1985) and Smith et al. 
(1988). Note that SNR and CR dynamics are not sensitive 
to the precise value used for T , because the shock struc- 
ture is mainly determined by the Alfvenic Mach number. 

We use an upstream magnetic field value B = 40 ^G, 
which is required to provide the required shape of the 
synchrotron spectrum in the radio and X-ray bands (see 
below) . 

The gas dynamic problem is characterized by the fol- 
lowing length, time, and velocity scales: 



R = (3M cj /4^ ) 1/3 , to = Ro/Vo, V = ^2E sn /M ei , 

which are the sweep- up radius, sweep- up time and mean 
ejecta speed respectively. Here po — lAm p Nn is the ISM 
mass density, m p is the proton mass. 

According to Chevalier ( 1982 ), an analytical approxi- 
mation to the shock expansion law during the free expan- 
sion phase (t < t ) is then 



R s o,Eg^/ 2k P - 1/k t^ 3 y^- 2 \ 

which for k — 7 gives 

R s oc {ElJ Po flh^. 

In the adiabatic phase (t > to) we have 

R s oc (£U/p ) 1/5 t 2/5 . 



(3) 



(4) 



(5) 



(2002) 



The obser ved mean expansion law of Tycho's SNR (Tan & 
Gull |1985|) is R s oc i" with p = 0.46 ± 0.02. This is consis- 
tent with the average over the oute r rim of the remnant, 
as obtained by Reynoso et al. ( 1997 ), even though the lat- 
ter authors found significant azimuthal variations of the 
expansion rate, indicating a substantial nonuniformity of 
the SNR environment, as mentioned above. On the whole 
therefore Ty cho's SNR should be ne ar the adiabatic phas e 
(Strom et al. |1982| ; Tan & Gull |l985j ; Reynoso et al. |1997| ). 

The calculations together with the azimuthally aver- 
aged experimental data are shown in Fig.^. An explosion 
energy E sn = 0.27 x 10 51 erg is taken (in addition to the 
above ISM density, and the ejecta mass) to fit the ob- 
served SNR size R s and its expansion rate V s . This value 
is again in basic agreement with the earlier determinations 
referred to above. 

According to Fig.|l|a Tycho is indeed nearing the adia- 
batic phase. To fit the observed radio spectral index (see 
below) we assume a proton injection rate 77 = 3 x 10~ 4 . 
This leads to a significant nonlinear modification of the 
shock which at the current age of t = 428 yrs has a total 
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Fig. 1. (a) Shock radius R s and shock speed V s ; (b) to- 
tal shock (er) and subshock (c s ) compression ratios; (c) 
ejecta (-E j), CR (£c) 5 gas thermal (E gt ), and gas ki- 
netic (-Egk) energies as a function of time. Scale values 
are R = 2.72 pc, V = 4402 km/s, t = 605 years. The 
dotted vertical line marks the current epo ch. T he observed 
size and speed of the shock (Tan & Gull 1985 ) are shown 
as well. 



compression ratio a = 5.7 and a subshock compression 
ratio a s = 3.5 (Fig.ft)). 

We should at this point ask the question, how sensi- 
tive these results are to the specific choice of the over- 
all parameters. These parameters are partly the result 
of analyses of the thermal X-ray emission. The corre- 
sponding X-ray models need to take into account of non- 



equilibrium electron temperatures and non-equilibrium 
ionisation states behind the outer shock as well as non- 
solar chemical composition in the ejecta which are in part 
heated by the reverse shock. As a consequence, different 
workers do not obtain identical results. The differences to 
our set of parameters is most pronounced in the results 



of Hamilton et al. (1986) who derived the different values 
E sn = 0.7 x 10 51 erg, iV H = 0.28 cur 3 , and d = 3.1 kpc. 
One way of comparing the consequences of such different 
inputs into our theory is to check on their internal consis- 
tency. Taking into account that the linear size R s oc d, 
then it follows from Eq. (4) that E sn oc N^d 1 ' 2 for 
a given source. Changing our parameters iVn and d to 
TVh = 0.28 cm~ 3 and d = 3.1 kpc, respectively, we obtain 
E sn — 0.57 x 10 51 erg for the Hamilton et al. case, not 
too different from their actual value for E sn . We can also 
discuss the change expected for the resulting 7r°-decay 7- 
ray integral flux . Since in the strongly nonlinear case, 
which is characterized by efficient CR acceleration, the 
number density of accelerated CRs at relativistic ener- 
gies is proportional to poV 2 (e.g. Berezhko et al. 1996 



Berezhko & Volk 1997| ), the occupied volume goes like 
i? 3 , and the 7r°-decay 7-ray production is proportional to 
the gas density po, we have F* oc N^R^V 2 d^ 2 . Inserting 
Hamilton et al.'s values for d and JVh, the 7r°-decay 7-ray 
flux would be decreased to 0.77 times our value given in 
Fig|j| below. Given our value of the magnetic field, we can 
do the same for the electron-proton ratio K cp , given the 
measured radio energy flux vS v (see Fig.||). The value of 
K cp would be lowered by a factor of 1.3 times our value. 
The IC 7-ray flux would not change at all for given ra- 
dio flux. The magnetic field value which we deduce has a 
different character. We therefore defer its discussion to a 
later point of this section. 

We now continue with the nonthermal aspects of 
Tycho's SNR. 

Given the proton injection rate the acceleration pro- 
cess is characterized by a high efficiency in spherical sym- 
metry: at the current time t/to — 0-7 about 20% of the ex- 
plosion energy have been already transferred to CRs, and 
the CR energy content E c continues to increase to a maxi- 
mum of about 43% in the later Sedov phase (Fig.lc), when 
particles start to leave the source. As usually predicted 
by the spherically symmetric model, such a CR acceler- 
ation efficiency is significantly higher than required for 
the average replenishment of the Galactic CRs by SNRs, 
corres pond ing to E c « 0.1E SD . As discussed above (Volk 
et al. |2002 , Berezhko et al. 2002 ), this requires a physi- 
cal renormalization of the number of hadronic CRs. We 
choose / re = 1/5 here, in agreement with the theoretical 
estimate above, and consistent with the average Galactic 
CR replenishment. 

With this renormalization the CRs inside Tycho's SNR 
contain at the current epoch the total energy 



E c » 0.2/ ro £ sn « 10 49 erg 



(6) 
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Fig. 2. The overall CR spectrum as function of momen- 
tum. Solid and dashed lines correspond to protons and 
electrons, respectively. 




Fig. 3. Spectral energy distribution of the electron syn- 
chrotron emission. The observed X-ray ( Allen et al. 1999 ) 
and radio emission (Reynolds & Ellison 1992) are shown. 



The volume-integrated (or overall) CR spectrum 



N(p,t) 



16ttV 



drr 2 f(r,p,t) 



(7) 



14 eV is the maximum CR energy 
is limited mainly by geometrical 



is given in FigJ|| (Note that only the nonthermal part of 
the spectrum is pictured. That is why both the proton 
and electron spectrum have a maximum above suprather- 
mal proton energies. As discussed above, the spectrum 
N(p) is of course smoothly connected with the thermal 
distribution, not given here.) Far above suprathermal en- 
ergies it has, for the case of protons, almost a pure power- 
law form N oc jT 1 over a wide momentum range from 
10~ 2 m p c up to the cutoff momentum p max = e max /c, 
where e max w 2 x 10 
(Fig.^). This value p, 
factors, which are the finite size and speed of the shock, its 
deceleration and the adiabatic cooling effect in the down- 
stream region (Berezhko 1996 ) . Due to the shock modifi- 
cation the power-law index slowly varies from 7 = 2.2 at 
P ^ m p c to 7 = 1.8 at p > 10 3 m p c. 

The shape of the overall electron spectrum N e (p) de- 
viates from that of the proton spectrum N(p) at high mo- 
menta p > p\ « 10 3 m p c, due to the synchrotron losses in 
the downstream region with magnetic field « 240 /iG. 
This field is assumed to be roughly uniform and equal 
to the potshock field B2: -Bd — B2 = ob-E>o, where 
cr| = 1/3 + (2/3)cr 2 for a strongly turbulent field. We 
shall approximately take <tb = cr in the sequel. 

The synchrotron losses become important for electron 
momenta greater than 



Pi 



1.3 



10" 



yr 



10 /iG 



(8) 



Substituting the SN age t = 428 yr into this expression, we 
have p\ w 800m p c, in good agreement with the numerical 
results (Fig.||). 

The shock constantly produces the electron spectrum 
f e oc p~ q , with 5 « 4, up to the maximum momentum 



Pmax which is much larger than p\. Therefore, within the 
momentum range p\ to p^ ax , the electron spectrum is 
/ e oc p~ 5 due to synchrotron losses, and the corresponding 
overall electron spectrum is N c oc p~ 3 . 

The maximum electron momentum can be estimated 
by equating the synchrotron loss time and the acceleration 
time, that gives 



Pn 



6.7 x 10 4 



10 3 km/s 



/ (*-l) 




/ cr(l + <T B <t) 


V Bo ) 



(9) 



At the current epoch V s ks 3100 km/s which leads to 
a maximum electron momentum p° lax ~ 6 x 10 m p c in 
agreement with the numerical results (Fig.||). Overall, the 
parameters K cp = 4 x 10~ 3 and B± = 0.5B d — 120 fiG 
provide good agreement between the calculated and the 
measured synchrotron emission in the radio to X-ray 
ranges (Fig.||). This concerns first the deviation of the 
radio-spectrum from the test particle form S v oc v~ ' 5 . 
At the same time, the steepening of electron spectrum 
at very high energies due to synchrotron losses naturally 
yields a fit to the X-ray data with their soft spectrum. 
Such a smooth spectral behavior is achieved in a 40 /iG 



upstream field (see also Berezhko et al. 2002). This is in 
clear contrast to the type of phenomenological fit by e.g. 



Reynolds ( |l998| ) for SN 1006. 

Since the total number of accelerated protons has to be 
renormalized by the factor / re = 0.2, the above number of 
accelerated electrons then corresponds to the renormalized 
parameter K cp = 2 x 10~ 2 . This turns out to correspond 
reasonably well to the canonical value of 10~ 2 of the elec- 
tron to proton ratio in the Galactic CRs. The total energy 
of accelerated electrons in Tycho's SNR is approximately 
E* = 2 x 10 47 erg, as estimated from the fit to the observed 
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Fig. 4. Observed total radio flux of Tycho's SNR 



(Reynolds & Ellison f 992) as a function of frequency, with 
model spectra superimposed. The solid curve represents 
our selfconsistent calculation, whereas the dashed curve 
illustrates the shape of the spectrum that would be ex- 
pected in the test particle limit. 



synchrotron emission (Fig|j) and our estimated value of 

In detail the radio data can be fitted with a power 
law spectrum S v cx v~ a , whose index a — 0.607 ± 0.007 



(Reynolds & Ellison 1992) is considerably larger than 0.5, 



as it would correspond to an electron spectrum N e cx p~ 
produced in the test particle limit by an unmodified shock 
with compression ratio a = 4. For our choice of the proton 
injection rate, r\ = 3 X 10~ 4 , the shock is significantly mod- 
ified by the backreaction of the accelerated protons (see 
Fig.[l]b) with a total present compression ratio a — 5.7 at 
the present epoch. At the same time low energy electrons, 
with momenta p < 10m p c (e c £ 10 GeV) which produce 
synchrotron emission at v < 10 GHz, are primarily ac- 
celerated at the subshock which only has a compression 
ratio a s — 3.5. Therefore these electrons have a steeper 
spectrum N c cx p~ 22 that leads to a radio spectrum 
S v cx v~ 6 that fits the experimental data quite well (see 
Fig.^, where the synchrotron flux at radio frequencies is 
presented). The fact that the observed value of the radio 
power law index a exceeds the value 0.5 is considered as 
an indication that the shock is significantly modified. Note 
that a = 0.6 is the average value of the power-law index 
within the frequency range shown in Fig.^. In fact due 
to the concave shape of the electron spectrum the index 
slightly decreases with increasing frequency, with a = 0.61 
and a — 0.56 at the lowest and largest frequency respec- 
tively. We emphasize that in our view not only a strongly 
modified shock, which produces a steep electron spectrum 
at energies e e < 1 GeV, but also a relatively high up- 
stream magnetic field strength Bq = 40 /iG, compared 
with typical ISM values Bq — 5 /iG, and a corresponding 
downstream value B± = 120 /iG are unavoidably required 
to have the energies of the radio emitting electrons in the 
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Fig. 5. Observed total integral hard X-ray flux of Tycho's 



SNR as a function of photon energy (Allen et al. 1999) 
together with the theoretical synchrotron spectrum. 



steep part of their spectrum e c < 1 GeV. According to 
the model calculations of Lucek & Bell ( 200C| ), the exist- 
ing ISM field can be significantly amplified near the shock 
by CR streaming to result in such values. An upper bound 
for the field strength is discussed below. 

We note that the necessity of strong nonlinear shock 
modification and large values of the magnetic field 
strength in young SNRs to reproduce their steep radio 
spectra was argued for the first time by Reynolds & Ellison 
(1992). The modified shock parameters and required mag- 
netic field £?d = 10~ 4 — 10~ 3 G which they derived for 
Tycho's SNR are consistent with our results. 

The large perpendicular downstream magnetic field 
B± = 120 /iG, which leads to a substantial steepening of 
the electron spectrum at high energies, also naturally pro- 
vides a smooth cutoff in the synchrotron spectrum S v (v) 
at frequencies v > 10 18 Hz, that correspond to X-ray pho- 
ton energies e v > 5 keV. As it can be seen from FigJ| 
the spectral shape of the calculated integral flux e v F v (e u ) 
is closely similar to what is observed (Allen et al. 1999| ); 
the X-ray data are just approaching the cutoff region. The 
small difference in the amplitude is not relevant, because 
it could be easily reduced by a fine tuning of parameter 
values which we do not attempt to achieve here. 

The calculated IC emission off the CMB nonther- 
mal bremsstrahlung (NB) and the derived 7r°-decay 7- 
ray fluxes are presented in Fig.^j, together with existing 
experimental data. Note that the number of accelerated 
protons, which produce the 7r°-decay 7-rays, is renormal- 
ized by the factor / re = 0.2 relative to the calculation in 
spherical symmetry. 

With this renormalization, the hadronic gamma-ray 
flux is just below the HEGRA upper limit in Fig. ||. The 
electron synchrotron - and thus also IC and NB - fluxes 
have to remain the same, fixed by the radio and X-ray 
observations. 
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Fig. 6. IC (dashed line), NB (dash-dotted line), and tt°- 
decay (solid line) 7-ray spectral energy distributions, as a 
function of 7-ray energy. The observed 3c 7-ray upper lim- 
its (W - Whipple (Buckley et al |199S| ), H-CT - HEGRA 
IACT-system (Aharonian et al. 2001)), and the 95% con- 
fidence HA - HEGRA AIROBICC upper limit (Prahl et 
al. 1997), are shown as well. 



According to the calculation, the hadronic 7-ray pro- 
duction exceeds the electron contribution by a factor of 
about 2 at energies e 7 < 1 TeV, and dominates at e 7 > 
10 TeV (Figjfj). The NB is always small above the GeV 
range. We cannot say much about lower energy NB be- 
cause the electron distribution is not well known at those 
lower energies, as discussed earlier. In the region of the 
7r°-decay bump around 67.5 MeV, the hadronic gamma- 
ray spectrum may be distiguishable from the electron IC 
spectrum. However, the corresponding low particle ener- 
gies are so far away from the "knee" in the Galactic CR 
spectrum around 10 15 eV that such a bump would not 
be a proof for Tycho's SNR being a typical source of the 
nuclear CRs. In addition, the steep background of the dif- 
fuse Galactic 7-ray flux would be hard to remove for an 
extended source, as most nearby SNRs would be. 

The 7-ray spectra produced by the electronic and 
hadronic CR components have rather similar shapes at 
higher energies 10 GeV ^ e 7 ?S 1 TeV due to the syn- 
chrotron losses of the electrons. Therefore, the only clear 
observational possibility to discriminate between the lep- 
tonic and hadronic contributions is to measure the 7- 
ray spectrum at energies significantly higher than 1 TeV, 
where these two spectra are predicted to be essentially 
different. The detection of a substantial flux at energies 
e 7 10 TeV would provide direct evidence for its hadronic 
origin. 

3.1. B-field amplification possibilities; upper bound 

Strong shocks like those of young SNRs have their diffu- 
sive properties dominated by self-excited waves from the 



proximation the mean square field fluctuation reaches a 
level < SB 2 > /4tt ~ (VJV s )l/2pV s 2 , where 1/2 P V S 2 is 
the upstream flow energy density in the shock frame, and 
V s /V a 3> 1 is the Alfvenic Mach number of the shock, with 
V a = B/jAirp) 1 / 2 (McKenzie & Volk |1982| , see also Drury 
P.983D - Nevertheless < SB 2 > /B 2 ~ V s /V a which implies 
that the rms wave field dominates the regular field. 

Linear Alfven waves constitute the long-wavelength 
elastic response of a magnetized plasma to a transverse 
perturbation. In reality it is not clear why the wave field 
could not grow even more strongly in the nonlinear regime 
due to a nonlinear evolution of the CR streaming in- 
stability. Such a "plastic" response (Volk & McKenzie 
1981) makes the distinction between a mean field and 



its fluctuations rather meaningless, effectively increasing 
the magnetic field in which particles undergo spatial dif- 



CR streaming instability (Bell 1978). In quasilinear ap 



fusion (Volk |1984p . Bell & Lucek Q2001| ) have constructed 
an interesting phenomenological model for such a nonlin- 
ear development, arguing for a strong reduction of the 
effective diffusion coefficient. For strong fluctuations the 
diffusion coefficient should come close to the Bohm limit 
K ~ l/3r 9 c, with r g = cp/eB, where B is the effective 
magnetic field strength. The maximum momentum to be 
reached in diffusive shock acceleration is then proportional 
to B for any given time the process is assumed to operate. 

How large the energy density of the effective field B can 
actually become is a complex question. However, an upper 
bound is given by imposing no other constraint than the 
requirement that the Alfven Mach number always remains 
large compared to unity, so that the shock continues to be 
strong and to produce the nonlinear field amplification 
in the first place. Requiring then V s /V a < 5 we obtain 
Bl h /8n = (l/25)l/2pV; 2 . Applying this to the shock of 
a type la SN, we take for the interstellar field B lSM — 
5p G, V s = 5 x 10 3 km/s during the sweep-up phase, and 
J B ISM /(4 7 rp) 1 / 2 ~ 14 km/s, using our value of the ISM 
density ambient to Tycho's SNR. Thus B uh /B ISM ~ 71 
as upper bound for B. 

Our estimate for the strength of B in this paper con- 
sists in comparing the observed synchrotron spectrum 
with a theoretical spectrum for the shock accelerated elec- 
trons. We conclude that the -B-value inferred is still well 
within this upper bound, especially considering that it is 
actually the mean square field value that counts. 

An interesting corollary concerns the upper bound p u h 
for the proton cutoff momentum which is achieved at the 
end of the sweep-up phase. Given that typical estimates 
yield p uh ~ 10 14 eV in SNRs for B = B ISM , we obtain for 
Tycho's SNR p u b — 7.1 x 10 15 eV which corresponds to 
the so-called knee in the Galactic CR spectrum. 

Even if we were to assume a much higher value 
B ism /(4tt / 9) 1 / 2 = 50 km/s for the ISM Alfven velocity, 
these conclusions would not change qualitatively. This 
shows that at least in type la SNe the knee energy is the 
upper bound for diffusive shock acceleration. 
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4. Summary 

The kinetic nonlinear model for CR acceleration in SNRs 
has been applied in detail to Tycho's SNR, in order to 
compare theoretical results with the recently found very 
low observational upper limit for the TeV 7-ray flux. We 
have used stellar ejecta parameters Af C j = 1.4M Q , k = 7, 
distance d = 2.3 kpc, and ISM number density Nh = 
0.5 cm -3 . A total hydrodynamic explosion energy E sn — 
0.27 x 10 51 erg was derived to fit the observed size R s 
and expansion speed V s which are determined by the ratio 
E^ n /Nn. Even though the distance to the object is not 
very well known, the set of parameters has been shown 
to be internally consistent, and the predictions for the 
radio the and 7-ray fluxes are quite robust with respect 
to different parameter values in the literature. 

A rather high downstream magnetic field strength 
£?d ~ 240 /iG and a proton injection rate r\ — 3 x 10~ 4 
are needed to reproduce the observed steep and concave 
radio spectrum and to ensure a smooth cutoff of the syn- 
chrotron emission in the X-ray region. We cannot exclude 
that the required magnetic field strength, that is signifi- 
cantly higher than the rms ISM value 5 /xG, might have 
to be attributed in part to its nonlinear amplification near 
the SN shock by CR acceleration itself. The evidence for 
efficient nucleonic CR production that comes from the ra- 
dio and X-ray data and leads to a strong shock modifica- 
tion, is even more definite for Tycho's SNR than in the 



case of SN 1006 (Berezhko et al. |2002[ ). 

We find that, after adjustment of the predictions of 
the nonlinear spherically-symmetric model by a physi- 
cally necessary renormalization of the number of acceler- 
ated CR nuclei to take account of the quasi-perpendicular 
shock directions in a SNR, quite a reasonable consistency 
with most of the observational data can be achieved. The 
resulting nonthermal electron to proton ratio turns out to 
be consistent with the observed ratio in interstellar space. 
The total 7-ray flux at 1 TeV (with the 7r°-decay compo- 
nent exceeding the IC component) comes out to be slightly 
larger than the most restrictive observational upper limit 
from the HEGRA experiment. Given the remaining un- 
certainties in the basic astronomical parameters of Tycho, 
in particular those coming from analysis of the thermal 
X-ray data, we do not consider this as a problem for the 
theory. It rather leads us to the prediction that detec- 
tors with several times higher sensitivity, like MAGIC or 
VERITAS in the Northern Hemisphere, should indeed de- 
tect this source above 100 GeV in 7-rays. 

The expected 7r°-decay 7-ray flux cx e" 1 extends up 
to <; 30 TeV, whereas the IC 7-ray flux has a cutoff above 
a few TeV. Therefore the detection of 7-ray emission at 
~ 10 TeV would in addition imply clear evidence for a 
hadronic origin. 
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